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By means of exact diagonalization we study the low-energy states of seven electrons in the lowest 
Landau level which are confined by a cylindric external potential modelling the rest of a macroscopic 
system and thus controlling the filling factor v. Wigner crystal is found to be the ground state for 
filling factors between v = 1/3 and v = 1/5 provided electrons interact via the bare Coulomb 
potential. Even at v = 1/5 the solid state has lower energy than the Laughlin's one, although the 
two energies are rather close. We also discuss the role of pseudopotential parameters in the lowest 
Landau level and demonstrate that the earlier reported gapless state, appearing when the short- 
range part of the interaction is suppressed, has nothing in common with the Wigner crystalization 
in pure Coulomb case. 
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I. INTRODUCTION 

After the Laughlin states (LS) were proposed as new ground states of strongly correlated 2D-electron liquid in 
external magnetic field they were intensively compared to the known ground states (GS), in particular with the 
Wigner crystal (WC), to understand the conditions and limitations of the experimental observation of the fractional 
quantum Hall effect (FQHE). In fact, even LS themselves may be called "liquid" only for sufficiently large filling 
factors v = 1/m (for m < m c sw 71 ||), as follows from the formal analogy between the LS and the two-dimensional 
one-component plasma at dimensionless temperature T — l/2m. At very small T the equivalent plasma undergoes 
Kosterlits-Thouless transition to the state with the finite shear modulus and should be rather viewed as a solid @. 

For the Coulomb system, however, the critical value, m c , is of academic importance only, since it is easy to prove 
that in the solid phase the Lauglin state differs qualitatively from the genuine GS. Indeed, in a solid with non- 
zero shear modulus and Coulomb interaction between the particles the sound dispersion law in magnetic field is 
tvk ~ fc 3 / 2 Q. Calculating the mean-square displacement in such a solid at T = one finds a convergent answer 
< (u(0) — u(R — ► oo)) 2 >^ J d 2 k/uk — * const. On another hand, the equivalence between the LS and the finite- 
temperature 2D plasma implies divergency of this correlator. [In solid plasma this divergency is due to the transverse- 
sound dispersion law: < (u(0) — u(R — > oo)) 2 >^ T J d 2 k/u\ ~ Tln{R)\. Thus solid LS maintains the topological 
order only, which is typical for the 2D solid with short-range interactions in magnetic field when ujk ~ fc 2 , in agreement 
with the fact that LS is a perfect trial function for GS of the system with short-range interactions. 

Having established that in the Coulomb system LS provides incorrect GS for large m > m c , one may further suspect 
that it may give the way to the Wigner crystal at much smaller m. Early variational calculations performed for the 
electrons at the lowest Landau level and neglecting Landau-level-mixing effects || gave strong evidence that WC has 
lower energy already at v = 1/7, which explained why there is no Hall conductivity quantisation at this filling factor. 
Mixing effects, which are very important in real systems, were taken into consideration in Refs. ||. It was found 
that virtual transitions between the Landau levels promote WC states and make LS unstable even at v — 1/3 for 
sufficiently large mixing parameter A = v 1 / 2 {e 2 /elff)/u) c , where Ih = (1/eB) 1 / 2 is the magnetic length in the external 
field B, e is the dielectric constant, and u> c = eB/m* is the cyclotron frequency for electrons with the effective mass 
m*. (We use units % = c = 1). For the A = case the results of Refs. ||,|| predict LS to be the ground state for filling 
factors v = 1/3 and 1/5. 

An essential drawback of previous calculations is their variational character. Since the difference in energy (per 
particle) between the LS and WC is only a few percent in Coulomb units e 2 /e/#, only filling factors v = 1/m with 
odd m, were the liquid GS is known reasonably well, were analyzed in detail. We simply do not know other liquid 
states with necessary accuracy to compare them with the variational functions we have at hand for WC. Moreover, it 
is impossible to use LS energies i?is(l/m) to derive Er,i qu id{v) for other filling factors by extrapolation, because of 
the cusps that must occur at filling factors where FQHE exists. 
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It is believed that due to this cusping down at simple rational v there may arise reentrant WC-LS-WC behaviour 
around v — 1/5 (and possibly v = 1/3). Available experimental data seem to give strong evidence that WC exists 
at v — 0.21 However, this conclusion is based entirely on the divergent, activation-type resistivity p xx — > oo 

at low temperatures, thus one may not rule out the possibility of explaining the data by impurity induced electron 
localisation. 

The other way to study this problem is by means of exact diagonalization. We are not aware of any systematic 
attempt to look for the WC states in numerical simulations of FQHE, although in [|ll| WC was suggested as a 
possibility to explain perculiar degeneracies in the numerical spectra as a function of inter-particle pseudo-potential. 
As we show below, the collapse of the LS to the gapless state at v = 1/3 found in |ll|] , has nothing to do with the 
transition to WC state. The best numerical calculations were done on a sphere for as many as N — 12 particles on 
25 orbits (l^] . The spherical geometry being perfect for the study of FQHE itself, is practically ineligible for looking 
at Wigner crystallization. There is little doubt that the optimal electron configuration in the WC is the triangular 
lattice, which is topologically prohibited on a sphere. 

To give further theoretical support to the idea that WC state may be GS of the 2D electron gas at v < 1/5, we 
calculated numerically GS and low-laying excited states for the cluster of 7 electrons in the lowest Landau level and in 
the confining potential. (The number of particles in the system must be 7,12 etc., depending on the sample geometry, 
to account for the hexagonal symmetry of WC). The confining potential was derived from the Coulomb interaction 
between an electron in the cluster and electrons outside the "first coordination sphere" (that is at a distance equal 
or larger then y/3a, where a is the atomic length in the triangular lattice). In fact, the radius, a, of this confining 
potential was our main variable determining an effective filling factor (a ~ v^ 1 / 2 ). It was found that WC state (see 
below the discussion of what has to be thought of as WC for the system of only 7 particles) is GS of the system 
for all filling factors between 0.4 and 1/5 except for the region 0.34 > v > 0.294 (including v = 1/3), where LS was 
essentially present in the structure of GS. These results strongly suggest (keeping reservations for possible finite-size 
corrections) that LS at v = 1/3 is very close in energy to WC, and that in a perfect Coulomb system at v — 1/5 the 
ground state may be WC. There is almost no doubt then that WC must exist between these two filling factors. 

Our calculation completely ignores impurities and Landau-level-mixing effects, as well as possible screening of the 
bare Coulomb potential, and the role of electron delocalisation in the direction perpendicular to the 2D plane. Some 
of these factors are very important in real systems, and, e.g., electron screening, may work in favour of LS. We believe 
that experimental data of Refs. ]7]-||, demonstrating FQHE at v = 1/5, can be accounted for along these lines. We 
are planning to study different pseudo-potentials in a separate paper. 



II. HAMILTONIAN AND NUMERICAL METHOD 



We study a 2D system of N e = 7 electrons in the magnetic field in a confining potential. Since we intend to 
model the macroscopic system in the finite-cluster calculation, we derive this confining potential as resulting from the 
Coulomb interaction between an electron in the cluster and other electrons outside the " first coordination sphere" . 
The unit length a{y) in the WC lattice is expressed through the electron density n e and filling factor v = 2irl 2 H n e as 

\/3 9 9 47T . . 

^-a 2 n e = loia 2 ^, (2.1) 
<Z ^V3 

(from now on we measure all distances in units of Ih, which is kept fixed in our calculation). To account for the short- 
range correlations we place surrounding electrons on the coordination spheres of the WC state. One may better view 
our system as originating from the classical crystal with atoms sitting on the m = orbits ^j(zj) — "i> m= o(zj — Zj), 
with the positions Zj forming an ideal triangular lattice corresponding to the filling factor v. We then allow for 
full quantum dynamics of 7 electrons with Zj = and | Zj \= a while keeping other electrons frozen, but coupled 
to the first 7 ones by Coulomb forces. Thus obtained confining potential is not spherically symmetric, and may 
cause transitions changing the momentum of an inner electron by Am = ±6n. Obviously, this coupling will promote 
"crystal" state for the central cluster. To avoid this shortcoming we ignore all these transitions, keeping only the 
diagonal part of the interaction; this procedure is equivalent to rotational averaging of the confining potential. 

We also ignore Landau level mixing, which means that Coulomb interaction is the only energy scale in the problem. 
To simplify the notation we measure all energies in units e 2 / (eln)- Working in the symmetric gauge A = l/2H(—y, x) 
we place the electrons on up to 31 orbits corresponding to the angular momentum states ranging from to = to 
m = 30. The maximum number of orbits in a given calculation was defined by the condition that occupation numbers 
< n(m) >, giving the probability to find an electron in the state with the orbital momentum to, stop changing (at 
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the level of 0.001), when the number of orbits is increased, and that the lagest momentum state be empty with the 
same accuracy. 

The starting Hamiltonian then can be written as 

H = ^ t Kni,m2,m3,m4 (! L 1 fl L ! < 1 m38m4 +^^1 "'o-m a m ; (2-2) 

tni ,7712 ,m>3 ,ni4 m 

where creates an electron in the state 

r m 

| m) = e-' z l / 4 . (2.3) 



In the symmetric gauge the interaction matrix elements in the first term in (2.2) satisfy the conservation law mi 
ni2 — m3 — ?Ti4 = and may be written as 

(_l)mi+m 3 <mi,m 4 <rn 2 ,rn 3 



v muma , mtim< - 4 ni/2 J2 J2 

A Uli=i " 4 w fc 1= o fc 2 =o 



x( _ 2) fe 1+ fe 2 C^C™*C£C™*I{MMM , (2.4) 

JO 27T 

where M = mi + TO2, Cj = — j)l), and Vc(g) gives the Fourier component of the pair potential. For the case 

of Coulomb interaction Vc(q) = 27r/g, and the final expression simplifies to 

I i\OTi+m 3 r= <mi,1Ti4 <m 2 ,m 3 
fmi,m2,m 3 ,m4 M 4 nl/2 2^ 2^ 

z KLLi=i"H-) kl=0 k2=0 

x( _ 4)fel+fe2 ^M-h-A&J-llHC^C^C^C^ • (2.6) 

To construct the confining potential one has to calculate the diagonal matrix elements for one electron staying on 
orbit | m) and the other electron staying on orbit ty m= o(z ~ Zj) (for the diagonal matrix element it does not matter 
whether the state ^f m =o( z — Zj) is defined in the same gauge as states | m) or obtained by gauge transforming the 
state | to = 0)). Thus we have 

VL MF) = £ V m ( Zj ), (2.7) 

3=2,3,... 

where the sum is over all coordination spheres starting from the second one, and 

V m (Z,) = (-ir2-™mlp J { ^ )[]2 fdq^- fl^^e""' J (« I *j I) ■ (2-8) 

Here Jo(q) is the Bessel function. In practice we constructed the confining potential by summing over all coordination 
spheres inside the radius 1001 h- 

Our diagonalization procedure is arranged as follows. For the calculation of the groundstate level we use the 
standard modified Lanczos method with the straightforward iteration procedure (see, for example, ]l2|), while for the 
calculation of the lowest excited levels we apply more sophisticated method ]13|] . The set of approximate eigenfunctions 



is reconstructed from Relay's tridiagonal matrix |14|, and the trial wavefunction is expanded in it. As is known, the 
set inevitably involves a substantial number of spurious states, due to numerical errors. These states, however, may 
be easily identified by their negligible contribution to the expansion of the trial wavefunction. Upon exclusion of the 
spurious states the set is subjected to the orthogonalization and correction by Newton method. The relative (with 
respect to a characteristic interlevel spacing) errors in the energy level calculation are typically of order 10~ 13 -j- 10 -11 
for the groundstate, and of order 10~ 9 -j- 10~ 5 for some ten first excited states. Since the Hamiltonian conserves 
the total angular momentum we take advantage of this symmetry to proceed separately for each angular momentum 
sector. 
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III. RESULTS FOR THE LOW-ENERGY STATES 



Before presenting our numerical results for the ground and first exited states in the cluster, let us first discuss how 
one may discriminate between the liquid and solid states in such a small system. The most obvious solution is to look 
at the pair correlation function. From the symmerty considerations we expect (at least for small filling factors) that 
one particle will be always staying near the origin, and the rest 6 particles will have their density distribution being 
peaked at a distance ~ a(v) apart. These particles are mutually correlated over the angle 9 between their coordinates 
on the "first coordination sphere". The appropriate pair correlation fuction thus can be defined as 

g a (6) = (a | ^{z 1 )^{z 2 )^{z 2 )^{z l ) \ a) , | z x | = | z 2 \= a , (3.1) 

where 9 — arg(zi) — arg(z2) varies in the interval (0,7r). In the solid state we expect three well-defined oscillations in 
g(9), while in the liquid these oscillations should be strongly damped. It is difficult to predict apriori the amplitudes 
of oscillations, but it is known (see, e.g., |^]) that pair correlations in the LS disappear very rapidly at v = 1/3 and 
v = 1/5. Our definition of g is not quite standard, but we believe that its qualitative behaviour is the same (we verify 
this point explicitly below). Anyway, the abrupt change of the ground state correlator gc{9) as a function of v is 
indicative of the solid-liquid transition. 

One may also expect some qualitative differences in the structure of the low-energy spectra of WC and LS. By 
construction, our Hamiltonian is cylindrically symmetric and conserves the total angular momentum M = ^2 m mn(m). 
In the solid phase of the macroscopic system this symmetry is spontaneously broken. For the triangular lattice under 
study the symmetry is broken by coupling momenta Mq ± 6n (where Mq is the ground state angular momentum 
and n is an integer). Thus in the solid phase we expect the states | Mq ± 6n > to form a subset of the lowest excited 
states well separated from the rest of the excitation spectrum in these sectors. There is no special reason to have the 
lowest excitations at Mq ± 6n in a liquid phase, nor should they have much lower energies than excited states with 

M = Mq. 

One remark is in order here. In a really macroscopic solid, the lowest states are those corresponding to the system 
rotation as a whole, with the energy going as E ~ {M — Mq) 2 /L 4 where L is the system's size. The crystal symmetry 
is not present in the structure of the spectrum explicitly, but it is important that the states, mixed by the symmetry 
breaking fields, are among the lowest ones. In the finite system of only seven particles we do not expect the spectrum 
to be quadratic in M — Mq, since this property in the rotating solid is achieved by creating extra zeros in the wave- 
function ^(zj) — ^(^i, z 2 , ■ ■ ■ , z/v) for fixed {zi, . . . , Zj—i, Zj+i, ■ ■ ■ , Zn}- This procedure may be too costy in energy 
in a small system. For N e = 7 rotation is equivalent to the correlated motion of six particles. The first rotating state 
which requires no extra zeros in ^(zj) is that with | M — Mq |= 6. 

The other point concerns the consistency of our procedure of controlling the filling factor according to equations 



(2.1) and (2.8). Since the confining potential is derived from the crystal state, a natural question arises of how 
good is this approximation for modelling a liquid environment. There is no doubt that at v = 1/3 the ground 
state is well described by LS with the ground-state angular momentum Mq — 3N e (N e — l)/2 = 63. No matter 
how trivial, this fact is not at all predetermined by the numerical procedure used. Its validity was confirmed in our 



calculations, thus demonstrating consistency between Eq. (2.1) and the effective filling factor. Similarly, we observed 



that Mq = 5N e (N e — l)/2 = 105 when v = 0.198 in Eq. (2.1). The consistency of our "mean-field" procedure follows 



also from the fact that for all v < 1/2 the position of the maximum in the particle density p(R) coinsides with a(v). 

In Table I we present our data for the ground state angular momentum as a function of filling factor. For v > 0.705 
the system is described by the IQHE state with occupation numbers rii — 1 for i = 0,1,..., 6. After drastic 
transformations in the range of densities between 0.705 and 0.46 GS evolvs into the state with well-defined pair 
correlation function g(9). We note, that starting from rather high filling factor 0.587 the angular momentum of the 
ground state changes by 6. Also, the lowest exited state is always in the sector Mq ± 6. 

To identify the nature of GS we present in Fig.l the plots of g{9) calculated for critical filling factors vm where 
Mq jumps. While going from Mq = 45 to Mq — 51 — > 57 the pair correlation function develops more pronounced 
oscillations. We naturally consider this evolution as formation of more rigid solid state order in the system, although 
the filling factor seems to be too large here to expect WC state in a macroscopic system. If we ignore for the 
moment what is happening in Mq — 63 then the "crystal set" may be smoothly continued to higher momentum 
states 57 — > 69 — > 75 — > 81 ... — > 111 resulting finally in a quite impressive "long-range" order (see Fig. 2). With all 
the reservations concerning small system size we have to conclude that WC has lower energy than LS in the range of 
filling factors between 1/3 and 1/5. 

We also observe a well-defined structure of "satellite states" | Mq ± 6n > in the energy spectrum for small v : for 
example, when Mq — 81 we find that £^75 — Eq and E^ — Eq are some five times smaller than the energy of the first 
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excited state in sectors M — 75, 81, 87 (see Fig. 3). Note also the remarkable similarity between the low-energy specra 
in the basic set ol states with M = Mq ± 6n. 

Clearly, the state with Mq = 63 is special in that its pair correlation function is more "liquid-like" than g{9) for 
both Mq = 57 and Mg = 69. As mentioned above LS at v — 1/3 has M = 63, thus irregular behaviour of the pair 
correlation function in this sector may be due to the change of GS from solid to liquid. This suggestion seems to be 
correct, because the calculated projection of the exact GS for v = 0.32, i.e., in the middle of the stability interval of 
the sector M = 63 (see Table I), on the Laughlin state is as large as < ^^g^ | >~ 0.934, and the ground state 

energy is extremely well approximated by the variational value E^g 3 ^ =< ^^g^ | H | >• In Coulomb units 

we find E^g 3 ^ — Eq 3 ^ — 0.0134, while the energy of the first excited state in the sector Mg = 63 is almost five times 
higher, E[ 6 3) - E { £ 3) = 0.0621. Furthermore, there is no pronounced satellite structure in the low-energy spectrum 
when Mq = 63. Surprisingly enough, the ground state wave function and gc(9) are rather different from | ^^g 3 ^ > 
and ghs{Q)- It is clearly seen in Fig. 4 that gLs{@) is almost flat for large 9 and shows no sign of pair correlations 
across the diameter of our system. These correlations are present in GS. Also, in Fig. 5 we plot the average occupation 
numbers < n(m) >, calculated in GS and in LS. We see that < n(m) > in LS has much smaller amplitude at m = 
and more shallow minimum. As one might expect beforhead, the central particle is not at all localized in the liquid 
phase. 

To clarify the nature of such differences, we construct another variational state, which may be regarded as solid, 
| § (63) >. Consider two nearest solid states, e.g., | *[J 5) > and | * (81) > at some 0.255 < v < 0.276. We notice that 
their distribution functions < n(m) > are very close in shape (see Fig. 6), with one particle being localised on orbits 
with small m (actually m = 0, 1; the sum of < n(m) > before the minimum is almost 1), and the other six particles 
occupying extended states with large m. When going from Mg = 75 to M — 81 the value of M = ^2 m m < n(m) > 
changes by 6 almost entirely due to the change of the occupation numbers of six particles on the first coordination 
sphere, i.e., < n/ 81 )(m + 1) >«< nP 5 \m) > for large m. Considering | ^( 81 ) > as rotating state with all the pair 
correlations being preserved, we may construct the variational state | <J/( 75 ) > close to exact | ^g 5 ^ > according to 
the rule 

| *< 75 ) >~ ]T Cf^al^al^ . . . aL 2 _i< I > , (m i+1 > rm) , (3.2) 

{mi} 

where J2i m i = 81, and Cj^j are the corresponding exact amplitudes of the expansion | vp( 81 ) >= 

S{ m ■} ^{m'} Ili a mi I >. Notice that the first particle keeps its states. This procedure is well justified because 
the first particle is separated from the others by a deep minimum in the distribution function (with < n(m) > close 
to zero in minimum, see Fig. 6). To estimate the accuracy of this procedure we project thus obtained variational 
state on exact | ^q 5 ^ > and find the overlap to be 0.995. We apply now this method to construct | <J>( 63 ) > from 
| \]/( 69 ) > obtained at v — 0.32, to obtain a solid-state trial wave- function. In Fig. 4 and Fig. 5 we show the pair 
correlation function and < n(m) > of this state. Finally, the solid-state variational energy turns out to be as good as 
E S 3) - E G 3) = °- 0075 and thc overlap with GS is < *( 63 ) | *g 3) >= 0.953 (even better than that of the Laughlin 
state!). 

From these data we have to conclude that the genuine GS in the range of filling factors 0.294 < v < 0.340 is a 
mixture of solid and liquid phases with comparable amplitudes. Not only these two quite different states strongly 
overlap with the ground state and almost minimize the energy, but also < 5>( 63 ) | ^^g 3 ^ >— 0.817. That large overlap 
is, of course, the finite-size effect. Obviously, under these conditions no definite conclusion about the true GS of the 
macroscopic system is possible, and there is no contradiction with the experimental fact that at v — 1/3 thc GS is 
the Laughlin liquid. 

We would like to comment here on the widly used argument, based on diagonalization of finite-size systems, 
that large overlap of LS with the exact GS and its precise energy, may serve as a criterion that the corresponding 
macroscopic system will be an incompressible liquid. We have demonstrated above that this argument simply does 
not work for the system of seven particles; short-range order in LS and WC turns out to be very similar. One has 
to analyze more delicate properties (like pair correlation function at large distances) to discriminate between the two 
phases. 

It follows from our data in Fig. 2 that GS in the sector M = 105 is of solid type. To see how different is | v[/^ 05 ^ > from 
I ^ls^ > we P resen t in Fig. 7 the corresponding correlation functions. We further confirm this result by calculating 
the overlap between the two states, < *^ 5) | *^ 105) >= 0.759, and thc Laughlin state energy E^ 5) - E { q 05) = 0.0188 
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(compare with the energy of the first excited state E[ ' ' — E G ' = 0.03f9). Now, the admixture of the LS in the 
structure of GS is much smaller than that at v = 1/3 and the variational energy is of the order of the first excited 
state in this sector. To reconcile this result with the experimental observation of the FQHE at v = 1/5 in some 
(not all!) systems [@-f§, we notice that our result was obtained on finite-size system and for the unscreened Coulomb 
interaction between the particles. Given rather large difference in energy between LS and WC found in our study, 
it is likely that WC will be the true GS of a macroscopic system too. This conclusion, however, may change for the 
screened Coulomb interaction since the Laughlin state is stabilized by short-range interactions. We plan to investigate 
the role of screening effects on the ground state at v — 1/5 in a separate paper. 

Since the liquid energy is casping down at v — 1/5, our results give very strong support to the idea that WC exists 
in the Coulomb system for v > 1/5. Even if WC is replaced with LS at v = 1/5 when the interaction potential is 
screened, it will most likely survive at slightly larger filling factors. We thus conclude that experiments |@-^ did see 
WC state around v — 1/5. 



IV. OTHER GROUND STATES IN THE PSEUDO-POTENTIAL APPROACH 

It was found in Ref. ]Tl|] that varying pair potential between the particles in the lowest Landau level one can 
drastically change the nature of the ground state. In this section we discuss whether this change is of any relevance 
to Wigner crystallization. 

Following Ref. |ll| we characterize the potential by the energies, U m , of pairs of particles with relative angular 
momentum m. In the lowest Landau level 

U m = jH dqq (YM^j e ~i 2 L m (q 2 ) , (4.1) 

where L m are the Laguerre polynomials. These are pseudopotential parameters because different bare interactions 
may have the same values of U m . For the Coulomb interaction these parameters are U m — y / 7r(2m — l)!!/(2 m+1 m!) 
and decrease slowly with m. Spinless fermions are coupled with odd values of m only. The effect of decreasing U\ 
for the Coulomb system of N e — 6 electrons on a sphere at v = 1/3 was the collapse of the Laughlin-type ground 
state to some gapless state [jllj (we will call it t^i-state). The nature of this state was not clearly identified, although 
the results did suggest a tendency to charge density wave formation. As we demonstrate below, the gapless ground 
state obtained by reducing the short-range part of the Coulomb interaction is not the conventional Wigner crystal 
(by "conventional" we mean the single- atom triangular lattice). 
We start by noting that the new state has almost zero overlap with LS |TlJ. This result is in sharp contrast with 



the large overlap between WC and LS found in Sec. [II. This fact alone is sufficient to rule out WC as a candidate for 
the C/i-state. Furthermore, as is seen from the data presented in Ref. p| , the collapse of LS is not accompanied by 
formation of the low-energy satellite states corresponding to the rotations of the octahedron formed by six particles 
on a sphere. 

We performed an analogous study of the ground state changes as a function of the U\ pseudopotential parameter 
for our system of seven particles. In agreement with Ref. Jll| we observe a drastic transformation of the ground state 
at v — 0.32 when U\ is reduced to 0.35. For smaller values of U\ the ground state angular momentum changes from 
Mq = 63 to Mq = 56. The change of Mq by 7, not by 6, also proves that we are not dealing with the conventional 
WC. Finally, we followed the transformation of the solid ground state with Mq = 75 at v = 0.265 and observed its 
collapse to the same t/i-state for U\ < 0.32. These results leave no doubt that reducing the short-range part of the 
Coulomb potential promotes new ground state other than LS or WC. 

To have a better filling about real-space interaction potentials with reduced values of U\ we show in Fig. 8 the 
particular set of interaction potentials of the form 

v{r) = yz e -v8 Jo(r 2 /8) _ A e _ rV4 1 (42) 

where I a is the Bessel function. The first term gives the Coulomb interaction between the two unit charges at a 
distance r = \r\ — r2 1 apart, each being spread out with the Gaussian distribution (27r) -1 / 2 exp{— | z — | 2 /2}, 
and the second term further suppresses the short-range part of the first. The choice of V(r) in this form is kind of 
arbitrary. It is justified by the simplicity of its Fourier transform V(q)/2ir — (1/q — A) exp{— q 2 }. In a more general 



case one may also vary the "cutoff length" by letting r — > r/r c in the second term. In Fig. 8 we plot the potential (4.2) 
for A = 0.8, 1.0, and 1.2. The corresponding values of U m are given in Table II. We see that C/i-state is stabilized at 
the edge of digging a potential well at short distances. 



6 



In Fig. 9 we present < n(m) > distribution in the t/i-state with momentum Mq — 56. Quite unexpectedly in the 
fJi-state, the central particle is replaced with the correlated hole. One has to appreciate this result in the system with 
the long-ranged Coulomb potential - by taking the central particle from orbits with m — 0, 1 and placing it to much 
higher orbits we substantially increase its mean-field energy. On another hand the " first coordination sphere" of six 
particles moves to internal orbits thus gaining some mean-field energy. Thus we see, that l^-state suggests locally ( 
on the scale of a) inhomogeneous particle distribution. Of course, the long-range tail of the Coulomb potential ensures 
that the macroscopic system is homogeneous on a large scale 3> a, but when the short-range part of the interaction 
is reduced, the system may choose states with local density higher than average. We are not able to say anything 
definite about such a state except that it is not conventional WC. Obviously, if the final state is a solid with more 
than one particle in the unit cell, it can not be traced from the numeric study of seven particles. 
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VI. NOTE ADDED 

After this work was completed we became aware of the fact that we had overlooked some important experimental 
results [ p"5| which seem to be in an excellent agreement with our numerical study. In these Refs. a metal - insulator 
transition is found to occur at the universal filling factor v c ~ 0.28 in rather wide range of magnetic fields and sample 
mobilities, no reentrant behavior is observed around v = 1/5. The authors argue that their results could be explained 
in terms of Wigner crystallization [though other interpretations are not ruled out] . 

TABLE I. GROUND-STATE ANGULAR MOMENTUM 

GS angular 
momentum 

M G 21 28 33 39 45 51 57 63 69 75 81 87 93 99 105 111 



1 0.705 0.587 0.527 0.460 0.408 0.364 0.340 0.294 0.276 0.255 0.240 0.224 0.211 0.198 0.188 
0.705 0.587 0.527 0.460 0.408 0.364 0.340 0.294 0.276 0.255 0.240 0.224 0.211 0.198 0.188 



TABLE II. PSEUDOPOTENTIAL PARAMETERS 



Potential Ui U 3 U 5 

Coulomb 044 028 022 

A = 0.8 0.37 0.29 0.23 

A = 1.0 0.35 0.28 0.23 

A = 1.2 0.32 0.27 0.22 
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FIGURE 1. 

Pair correlation functions g{0) for the ground states in the degeneracy points corresponding to the angular momentum 
changes 51 -> 57, 57 -> 63, 63 -> 69, and 69 -> 75 

FIGURE 2. 

Pair correlation functions for the ground states at the degeneracy points corresponding to the angular momentum 
changes 87 -> 93, 93 -> 99, 99 -> 105, and 105 -> 111 

FIGURE 3. 

Low-energy spectrum a,t u = 0.248 (Mg = 81) 
FIGURE 4. 

Pair correlation functions g(6) &t v = 0.32 for the ground state, the Laughlin state | \&£g 3 ^ >> an< ^ tnc sou d state 

| vp(63) > 

FIGURE 5. 

< n(m) > distributions at v = 0.32forthegroundstate, the Laughlin state | ^^s^ >> anc ^ the son d state | \]>( 63 ) > 
FIGURE 6. 

< n(m) > distributions for the ground state | ^q 5 ^ > and the excited state | <3/( 81 ) > a t ^ = 0.265 
FIGURE 7. 

Pair correlation functions g{6) at v = 0.193 for the ground state and the Laughlin state | ^ls^ > 
FIGURE 8. 

Some realizations of the interaction potential in the real space with reduced values of U±. 
FIGURE 9. 

< n(m) > distribution for the t/i-state at v = 0.32 and A = 1.2. 
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